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Ramona Gr6ber/’[^ Margarete Miihlleitner,^’|^ Eva Popenda,^’]^ and Alexander Wlotzka^’j^ 

^INFN, Sezione di Roma Tre, Via della Vasca Navale 84, I-OOI 46 Roma, Raly 
^Institute for Theoretical Physics, Karlsruhe Institute of Technology, 

Wolfgang-Gaede Str. 1, D-76128 Karlsruhe, Germany 
^ Paul Scherrer Institute, GH-5323 Villigen PSI, Switzerland 

We investigate the decays of the light stop in scenarios with the lightest neutralino being the 
lightest supersymmetric particle, including flavour-violating (FV) effects. We analyse the region 
where the three-body decay h —>■ Wbxi is kinematically allowed and provide a proper description 
of the transition region between the three-body decay and the four-body decay ii —> Xibff'. The 
improved treatment has been implemented in the Fortran package SUSY-HIT and is used for the 
analysis of this region. A scan over the parameter range including all relevant experimental con¬ 
straints reveals that the FV two-body decay into charm and x? can be as important as the three-, 
respectively, four-body decays if not dominant and therefore should be taken into account in order 
to complete the experimental searches for the light stop. 


I. INTRODUCTION 


The discovery of a new scalar particle by the LHC ex¬ 
periments ATLAS [T] and CMS [5] has marked a mile¬ 
stone for particle physics. The immediate investigation 
of its properties allowed to identify it as the Higgs bo¬ 
son, i.e. the quantum fluctuation associated with the 
Higgs mechanism. But still, the question remains open 
if it is the Higgs boson of the Standard Model (SM) 
or of some new physics (NP) extension beyond the SM 
(BSM). Among the numerous NP models that are inves¬ 
tigated, supersymmetric theories mn] certainly belong 
to the best motivated and most intensely studied BSM 
scenarios. Based on a symmetry between fermionic and 
bosonic degrees of freedom each SM particle has a super- 
symmetric (SUSY) counterpart. The SUSY partners of 
the top quark, the stops, play a special role. The large 
top quark mass allows for a large splitting between the 
two stop mass eigenstates H and ^2 with interesting phe¬ 
nomenological consequences. Thus, while the limits on 
the squarks of the first two generations are pushed to 
higher and higher values imiini, light stops have not 
been excluded yet by the experiments. Stops play an 
important role in the corrections to the SM-like light 
Higgs boson mass of the Minimal Supersymmetric Ex¬ 
tension of the SM (MSSM) and are crucial to shift its 
mass value from the tree-level upper bound given by the 
Z boson mass to the experimentally measured value of 
^ 125 CeV. Naturalness arguments favour the stops to be 
light as they significantly drive the amount of fine-tuning 
at the electroweak scale |1D]. In the MSSM with five 
Higgs bosons, two neutral CP-even ones, h and H, one 
neutral CP-odd one. A, and two charged scalars H^, the 
maximal mixing scenario optimally reduces the amount 
of fine-tuning m while ensuring the correct mass value 


* groeber@roma3.infn.it 
t milada.muehlleitner@kit.edu 
t eva.popenda@psi.ch 
5 alexander.wlotzka@kit.edu 


of h. Furthermore, light stops can help for the correct 
relic density through co-annihilation in scenarios with 
small mass differences of 15-30 CeV between the light 
stop and the lightest neutralino Xi P21 - I57] . And last but 
not least, light stops are necessary for baryogenesis to 
generate the matter-antimatter asymmetry in the MSSM 

EBHIO]- 

There exist numerous experimental analyses searching 
for stops in different mass windows. Light stops with 
masses below the kinematic threshold for the decay into 
a top quark and the lightest neutralino, assumed to be 
the lightest SUSY particle (LSP), can decay through the 
three-body decay U —into the LSP, a W bo¬ 
son and a bottom quark b. If the U mass lies below the 
three-body decay threshold, the light stop, assumed to 
be the next-to-lightest SUSY particle (NLSP), can de¬ 
cay through a FV process into the LSP and a charm 
quark c or an up quark u, ti —)■ (m/c)xi |41I |42]. An¬ 
other competing decay channel in this mass regions is 
the four-body decay H —>■ Xibff |13], where / and /' 
stand for generic light fermions. Former bounds on the 
stop masses have been set by LEP |44l l45] and Tevatron 
111113 . Searches based on charm tagging and monojets 
have been performed by ATLAS |1S] and CMS [11 . More 
stringent bounds have been derived by ATLAS in decays 
into charm quarks or in compressed SUSY scenarios in 
m as well as in final states with one isolated lepton, jets 
and missing transverse momentum m- ATLAS searches 
in final states with two leptons have derived bounds on 
the stop mass under the assumption that it decays into 
a 6-quark and an on-shell chargino, which decays via a 
real or virtual W boson, or that the stop decays into a 
top quark and the lightest neutralino [ 52 ] . The same de¬ 
cay modes have been taken in the analysis performed by 
CMS [ 53 ]. The latter analysis provides limits for various 
assumptions on the branching ratios, while the former 
analyses assume branching ratios of one in the respective 
final states. 

In |54| we have reinterpreted the charm-tagged and 
monojet searches gMIl by taking into account that the 
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branching ratios for the FV two-body and for the four- 
body decay can deviate significantly from one. This leads 
to considerably weakened exclusion bounds. In this work 
we investigate the transition region at the threshold of 
the three-body decay ti —)■ Wbxi- In particular, we 
analyse in this threshold region the interplay between 
the FV two-body decay and the three-body decay above 
the threshold, respectively, the four-body decay just be¬ 
low the thresholdjj It turns out that the two-body decay 
can still be significant here for certain parameter config¬ 
urations and can hence be exploited to improve and/or 
complement analyses based on the three-body decay fi¬ 
nal states. We extend and refine former analyses [55H57] 
by including the recently computed SUSY-QCD correc¬ 
tions to the FV two-body decay |540 and the FV tree- 
level couplings in the three-body decay as well as in the 
four-body decay where also a non-vanishing r and bot¬ 
tom mass in the final state jSl] are taken into account. 
Furthermore, the transition region between three- and 
four-body decays is consistently described by including a 
finite width in the W boson propagator, which becomes 
virtual below the three-body decay threshold. Finally, we 
check for the accordance with the LHC data on the Higgs 
boson search, the exclusion limits from SUSY searches as 
well as constraints from the relic density and H-physics 
and from electroweak precision measurements. 

In Sec. El details on the calculation of the three- and 
four-body decay widths are given, followed in Sec. |HI| by 
the description of the numerical analysis and the applied 
constraints. Our results are presented in Sec. |TV] We 
conclude in Sec. (V] 


II. THREE- AND FOUR-BODY STOP DECAYS 


proceeds via down-type squark, chargino and up-type 
quark exchange. The Feynman diagrams are displayed 
in Fig. [T] The index s = 1,..., 6 of the exchanged squark 
refers to one of the six squark mass eigenstates, which 
are not flavour eigenstates any more. In case of small 
FV as given in the MFV setup, the dominant final state 
is given by Wbj('(. We have calculated the three-body 
decay with the general flavour structure by extending 
the results of [53 [57] to all flavours. The full depen¬ 
dence on the bottom quark mass has been taken into 
account, whereas the first and second generation quark 
masses have been set to zero. The result for the decay 
width has been checked against a second, independent 
calculation by using FeynArts/FormCalc [HH - iHT] . 

In the threshold region where the three-body decay 
mode of the light stop into Wbxi opens up, the off-shell 
effects of the W boson can be described by implementing 
the W boson width in the propagators of the W boson 
diagrams in the four-body decays 

ui ^ X?d./f ■ (2) 

Again in case of small FV the dominant final state is 
the one involving the 6-quark, i.e. di = b. The W boson 
width in the propagators introduces a gauge dependence. 
The width renders the W boson mass mw in the W boson 
propagators complex, whereas it is real in the correspond¬ 
ing Goldstone boson couplings, so that the cancellation 
of the gauge parameter dependence between the W bo¬ 
son and the associated Goldstone boson diagrams cannot 
take place any more. Possible solutions are given by the 
complex mass scheme |5S] , where a complex mass is intro¬ 
duced also in the Feynman rules, or by the overall-factor 
scheme [ 511170 ], in which the whole tree-level amplitude 
is multiplied by 


We work in the framework of the MSSM with gen¬ 
eral flavour structure. Flavour-violating effects are con¬ 
strained by experiment to be very small which can be 
naturally accounted for in the Minimal Flavour Violation 
(MFV) |591l58] approach e.g., where the only sources of 
FV are given by the GKM matrix elements. Flavour vio¬ 
lation is induced through renormalisation group running. 
Due to the large mixing in the stop sector, the lightest up- 
type squark ui is then mostly stop-like. For convenience, 
we occasionally refer to it as the light stop U in the fol¬ 
lowing although it is understood that it has small flavour 
admixtures from the charm- and up-flavours. The three- 
body decay of ui into the lightest neutralino, a down-type 
fermion (* = 1, 2, 3), where i denotes the quark flavour, 
and a W boson. 


Ml Wdix\ , 


( 1 ) 


^ Note that we choose the parameters such that in the four-body 
decay only the diagrams with the intermediate W boson can 
become on-shell in the investigated region. 

^ See also m- 


n 

W propagators 


Pw - 

Pw ~ ’ 


(3) 


where pw denotes the W boson four momentum and T w 
the W boson width. The product ]/[ accounts for the 
maximal number of W propagators in the amplitude. 
We use the overall-factor scheme to ensure a gauge in¬ 
dependent result. The drawback of this method is that 
close to the threshold the non-resonant contributions are 
neglected. We checked, however, explicitly, that in the 
scenarios found in our numerical analysis below, the ef¬ 
fect of neglecting the non-resonant contribution is less 
than about 2% and hence acceptable. The three-body 
decay and the thus calculated four-body decay widths 
have been implemented in the SDECAY [73 CD routine 
of SUSY-HIT [73|, where the SUSY Les Houches Accord 
(SLHA) [71] format has been extended to the SLHA2 for¬ 
mat [75], as described in Ref. [M], to account for FV. 

In order to ensure that the three-body and the four- 
body decay widths match for mass differences 

above the kinematic threshold of an on-shell W boson, 
the W boson width must be computed in accordance with 
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the loop order and the input values used for the compu¬ 
tation of the four-body decay width. Thus, the tree-level 
W boson decay width is computed with massless first 
and second generation fermions, while the masses of the 
bottom quark and the r lepton are kept finite. 

III. NUMERICAL SETUP AND 
EXPERIMENTAL CONSTRAINTS 

We have performed a random scan over the parameter 
space of the model with the same settings as in the U ( 2 )- 
inspired scan of Ref. [M]. The parameters have been 
varied in the ranges 

tan j3 G [1,15] , 

Ma G [150,1000] GeV , 

Ml e [75,500] GeV , 

M^^ € [300,600] GeV , 

Mq^ G [1000,1500] GeV , 

At G [1000,2000] GeV . (4) 

The remaining MSSM input parameters have been chosen 
as 

M 2 = 650 GeV , 

M 3 = 1530 GeV , 

H = 900 GeV , 

Mp, =Mf =1 TeV , 

^1,2,3 -^1,2,3 ’ 

Mn = Mfj = Mp = 1.5 TeV , 

Wl,2 ^1,2 -^1,2,3 ’ 

Au = Ae = Ae = 0 . (5) 

The SM input parameters have been set to the PDG 
values m- We have applied the same constraints on the 
generated parameter points as in |54[ . but updated the 
branching ratio of the -G pT decay to the recently 
reported value 

^ m+M”) = (2.8;°;^) X 10-9 [77] . (6) 

Additionally we have checked for the dominant restric¬ 
tions due to electroweak precision observables by throw¬ 
ing away all points which are outside the 2a interval 


around the experimental value for the p-parameter 

p= 1.0004 ± 0.00024 IZg . (7) 

Among the parameter points fulfilling the constraints we 
have retained only those, for which the masses of the 
lightest up-type squark ui and the lightest neutralino Xi 
comply with 

— m^o € [60, 140] GeV . ( 8 ) 

The mass window around the three- to four-body de¬ 
cay threshold has been chosen large enough to allow for 
studying all effects that emerge in the threshold region. 
Finally, for the parameter points above the threshold 
SModelS |78l[Hn] based on the tools Phythia 6.4 |81| . 
NLL-fast |821155] and PySLHA [SS], is used to ensure that 
all parameter points fulfil the exclusion bounds derived 
from direct searches by ATLAS and GMS |49ti551 [50]. 
Since the searches in the FV two-body decay channel are 
not covered by SModelS yet, for the parameter points be¬ 
low the threshold the procedure explained in is used. 
The scenarios surviving all constraints include chargino 
masses around 660 GeV, slepton masses of 0(1 TeV) and 
charged Higgs masses in the range ^ 400 to ^ 1 TeV, 
so that the corresponding diagrams in the four-body de¬ 
cay with these particles in the intermediate propagators 
never go on-shell in the investigated threshold region. 

IV. RESULTS 

Figure shows the two-, three-, and four-body decay 
widths, respectively, for the parameter points of our scan 
which are in accordance with all applied constraints. The 
three-body decay, given by the green points, sets in at the 
the threshold rriu, — m^o = mw ± to^.. As expected, it 
approaches the four-body decay, illustrated by the blue 
points, for rriu^ — m^o mass differences sufficiently above 
the threshold^ The relative size of the four- and the 


^ Note, that we have implemented the total width of the top quark 
in the three-body decay and explicitly checked that the top width 
does not play a role in the three-body decay, also for rriui 
mass differences as large as 140 GeV. The three-body decay with 
the top quark width and FV couplings has been implemented in 
the SUSY-HIT version that includes FV decays. 
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FIG. 2. Partial widths of the ui two- (red), three- (green) 
and four-body (blue) decays as a function of the — m^o 
mass difference. 



FIG. 4. Branching ratio of the u\ two-body decay as a func¬ 
tion of the rtiux difference. The dashed line marks 

the threshold for the three-body decay. 


three-body decay widths is displayed in Fig. It shows 
that the finite width effects are still sizeable 20 GeV above 
the threshold and therefore should be taken into account, 
as is done by including the total width of the W boson 
in the four-body decay. Note, that the scattering of the 
points at the upper end of the mass difference is subject 
to the numerical integration precision in the four-body 
decay. Furthermore, the remaining off-set between the 
four- and three-body decay at large mass differences is 
due to the finite value of the W boson width. 

As can be inferred from Fig. the values of the two- 
body decay widths are equally distributed along the cho¬ 
sen range of the —m^o mass difference. While above 
the threshold the three-body decay dominates, close to 
the threshold the decay width for the two-body decay, 
which is shown in red, can be of similar size as the three- 
and four-body decay width, respectively. The branching 
ratio of the two-body decay is depicted in Fig. With 
possible values as large as ~ 40% at 20 GeV above the 


[-H 


[-H 



m^j — m^o [GeV] 


FIG. 3. Comparison of the ui four- and three-body decay 
widths as a function of the ma, — m^o mass difference. 


threshold, the two-body decay clearly is competitive with 
the other decay modes and thereby offers new discovery 
perspectives for light stops in this parameter region. In 
this region the charm is not soft any more and charm 
tagging could be used efficiently, as has been shown in 
m, where a search for pair produced scalar partners of 
charm quarks was performed. Such large two-body decay 
widths are achieved in scenarios with relatively large FV 
as is the case in the U (2)-inspired scenarios investigated 
here. If such a set-up is realised by nature, Fig. shows 
that it might not be possible to detect the light stop 
in the three- and four-body decay mode, respectively, if 
the masses of the light stop and the neutralino are such 
that they fall into the threshold regime. Hence, comple¬ 
mentary searches in the two- and the three-, respectively. 
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FIG. 5. Parameter points of the scan, surviving all applied 
constraints, in the “^xj plane. The colour code indicates 
the corresponding values of the FV two-body decay branching 
ratios. The upper (lower) grey line marks the lower (upper) 
bound of our investigated mass interval below (above) the 
threshold for the three-body decay. 
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four-body decay mode are required in this case. 

The two-body decay branching ratios for all scenarios 
of the random scan that passed the constraints are plot¬ 
ted in Fig.ISlin the rriu^-m^o mass plane. The upper and 
lower grey mes mark the borders of the interval defined 
in Eq. ([^ and the colour code indicates the value of the 
two-body decay branching ratio. While for the low mass 
region the parameter space is already very constrained 
such that no valid parameter points with stop masses 
lower than 260 GeV have been found, the fade out at 
high values of the stop and neutralino masses in Fig. is 
due to the limited scan range of the input parameters of 
the model. The plot nicely illustrates the relative impor¬ 
tance of the two-body decay in the four- to three-body 
transition region and underlines once more the necessity 
to take this decay channel into account in order to al¬ 
low for a proper analysis of the stop decays in this mass 
range. 


V. CONCLUSION 


Appendix A: Three-body decay width 

In this appendix we give the analytic formula for the 
decay width of the process ui —?> WdiXi- We first define 
the couplings relevant for the decay width. The coupling 
of an up-type and down-type squark to the W boson, 
with the corresponding part of the Lagrangian given by 
C = defined as 

9w^J. , (Al) 

^ bi=i 

where denotes the CKM matrix, kF“ and the 

squark mixing matrices in the SCKM basis as defined 
in the SLHA2 convention m and 52 the SU{2)l cou¬ 
pling constant. The squark-quark-neutralino couplings 
to up-type quarks and squarks, with C = + 

9xuiu,PR)'^sXii are defined as 


In this paper we have analysed decays of the light stop 
including the possibility of FV and with the lightest neu¬ 
tralino being the LSP. We investigated in particular the 
mass range where the three-body decay into Wbxi is 
kinematically allowed. We provide a proper description 
of this threshold region by resorting to the four-body de¬ 
cay into Xi&//' where the W boson total width has been 
taken into account in a gauge invariant way. The result¬ 
ing decay formula and the three-body decay with general 
flavour structure have been included in SUSY-HIT. We 
performed a scan over this threshold region where only 
the points in accordance with the constraints from the 
LHC Higgs and SUSY data, from the relic density and 
B physics measurements as well as from the electroweak 
precision data have been retained. The investigation of 
these scenarios revealed that the FV two-body decay into 
cx5 can be comparable to the three-, respectively, four- 
body decay and even dominate for some parameter sets. 
In order to properly investigate this mass region, the ex¬ 
periments should therefore also investigate two-body de¬ 
cays with charm quarks in the final state, in order not 
to miss the light stop, which might be the first SUSY 
particle to be discovered at the LHC. 
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'/2mw sin fj 
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Here, the W boson mass is denoted by mw, gi is the 
t/(l)y coupling constant and Z the neutralino mix¬ 
ing matrix as defined according to the conventions in 
Ref. . The angle /3 is determined through tan /3 = 
Vul^d where Vu^d denote the vacuum expectation values 
of Hu^d, respectively. The squark-quark-chargino cou¬ 
plings, with C = di + g^+^^-^PL^ UsXT^ are 
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given by 


92 


y/2mwsin, 

3 


i=i 


(A 6 ) 


-g2VnY.w:;vf, 


CKM* 


i=i 


9\,~ = 


92 


X, diUs ^/2mw cos j3 




i=i 


with U and V denoting the chargino mixing matrices as 
defined in El- The chargino-neutralino-lT boson cou¬ 
pling given by £ = xr(ffv+v Pl)XiW-, is 

defined as 


4u, = » 

Furthermore we introduce 


(A 8 ) 

(A9) 


(AlO) 


Note that we will drop the indices for Xi and denote the 
corresponding p. with By di we denote the final state 
down-type fermion with flavour index i. In SUSY-HIT we 
have set the masses of the first and second generation 
quarks to zero, corresponding to 2 = 0 and /i„i 2 = 0 . 
The differential decay width can then be written as 


dr = _”!.“\„ Re 


\Mf + \M^^\^ + \Mu\^+ 


32(27r)3 

2MjM*+ + 2M^Ml + 2M^+Ml 


dxidx2 


(All) 


and needs to be integrated over the reduced energies of 
the final state particles xi = 2{pux ' Pdi)l'rn\_^ and X 2 = 
2(Pfii ‘ Px)l'^lii Px denoting the four-momentum of 
particle X. The integration limits for the reduced energy 
xi are given by 


(^l)min — 2y^JId^ (A12) 

= l + (A13) 


For a given value of xi the range of X 2 is determined by 

(Ew + Ed^)^ 


(^2)niin — 1 “1“ f-^x 


ml. 


(A14) 


Ew-m 


W ’^di) 




(^2)max — 1 T Px 


{Ew + Ed,)^ 

ml. 


(A15) 


Ely-m 


2 

W ' 


El - 


TOt 


m~. 

Here Eyy and Ed, are the energies of the IF boson and 
the down-type quark di evaluated in the rest frame of the 
incoming stop and neutralino, 

^ ml,-ml^xi+ml-ml+mly ^ ^ ^ 
Ew = I —— (Aibj 

2jml^ -mlxi+ml^ 


Ed^ = 


m: 


\,Xi - 2m] 


2,/m? -mlxi+ml 


(A17) 

"Ml ' "'‘di 

The individual contributions to the differential decay 
width in Eq. (All) read 


® 9wuid^ 9 


W u\dt 


^9xdid,3^d,dt 


yi[Pwiy 2 + 2 / 3 )^ - Px- Pd, - 2yi] 


(l-x:i+ Pw - Pd,)^^ -X:i + pw - Pd,) 

+ (9^d,dJxld,+ 9^d,dAldjVJ^ 

22/1 +Px+ Pdj - Pwiy 2 + 2 / 3 )^ 

{I-X 3 +PW - -X 3 + PW - Pd,) 


(A18) 


\Mu\^=J2 


„2 -\/CKM-\/CKM^ 
y2 ^kn ^jn 


2(1 X2 + Px Puk){l X2 + Px pUj) 

I — 2 (yyjPuk9xukui9xl,ju, + ^ j)) 


V Px {pdi + ^ 2/2 + 2 .y 2 pw ) 

+ 2 9xnkU,9xljUiV9ukPu, (2/1 + ‘^y2y3Pw) 

F 2 9xuku,9xuju, iy^ipd, Pw F 4 ^ 2 ) F 2 y3pd, 

F 4 y22/3 F Pwi^yivl - 22/22/3/^dJ) | , 

(A19) 
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with = /i „3 + iy/JI^Ttlrriu^ and F* denoting the top 
width. For t = 1, 2 the {1^ are equal to Hm- 




^ (^ — 'r-i — 




( 


^Z/* 




g . g +,. g + g + 

•^Xfc diUi^xZXi xTXI 


+ (i? 0 L)^ 


4^3 (yi + y2 + + yi(Mx - t^w) 


+ 2AixJ/2(l - t^wVs) 


(g\, _ . - g\ g^X + (i? O L)) 

V x7di“i^xr^i“i xixi^xrxi / 

V^xt^^xt +2^^^ 2 / 22 / 3 ) 

( 2 ?^+. _ 27^: , - 2 ?^+ + (i? O L)) 

V XfcXi'^xrxi ^ 


y/J^+ik ^ 1) 

( 


3v^(yi + 2 / 2 ) 






g . g +,. g + g + 

•^Xfc Xfc Xi xr XI 


+ (i? -f-)- F)^ 


yJVXlJj^ + (/c O /) 


(Px + ‘^Vaf^w + 32 / 3 ) 


9^+h,,~,, 9-^+-^, 9^ + 


xtPui^xXdi'lii^xtxi^xTxi 4" ^ '^)) 

(^X^d.iil^Xi^d.iil^xtxi^Xi^Xl + ^ ^)) 

^y/JkdJPc if^W + Mx + 22 / 3 )I • 


(A20) 


The interference terms read 


6 2 


A4,A4*h 


-4 


dWuids 


=yy _ 

^ “ ^3 + - a;i - 

I (^x;^d.fil^x;^xi^xd.d, + ^ ■^)) 

y/Jhjy (2/1 - y2^^w{y2 + 2/3) + Md.) 

+ Xi4d.d. + ^ ^)) 

( 2/2 + 2/3)(Mx2/2 - 2 yiy 3 )g:^ + yi{2yi + 2/2 

- 2/3 + Mx) + Mx 2/2 - Md, (Mx + 2 / 3 ) 


2/3( 2/2 + 2 / 3 ) - Mx - yi) 


VJPkJPi ( 


„d?* „R* „L 

^xt diHl^xt Xl^xdid, 


+ {R -f-/- F)^ 


ted...3xf x.4.d, + ^ 4.) 

(A21) 


6 3 


-2y/2 32 


iti ds 


MrMi^yy -— 

jit ^ “ ^3 + Mx - /^dj(l - 3^2 + Mx - Ml. 


^CATM* ) /7^-^ R* _ L 
Xo^ y/Rujldx 9xujUi9y^did, 


• jn 


yi + Md, - y2 


yw (y2 + yz) 


+ 0^* - ~ 
w ax“3 «l^Xdid, 


3132 (^1 + 23 ^^^ 


(y2 + ys)) + yxy2 - yiy3 - 2.yl - yigd. 


+ y-d, (mx - y3) + yw yd. {-y 2 yz - yl) 


+ 9xuj fi 1 y id, d, x/y-<ii y" 


^X^jUl^Xd.ds 

ywy^iy^ + y^) 
yi +Md, -/^iv^y 2 (y 2 + y3) 


yi - Mx 


+ 9xUjUi9^did.,y^pk^ 


(A22) 


2 3 


52 


—1 v2(l Xi /d^+)(l X2 + /ix Mmj) 


yrCKM* I L R L* r—— - 

1 ^X^di«iyx^Xl^X">i V^^^xt 
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6nw + 6^3 


The notation (i? o L) means that the respective term is 
obtained from the previous one by replacing i? o L in 
the couplings, whereas {k y-)- 1) and {k ■(r^ j) means that 
the term is obtained by interchanging indices k and I or 
k and j, respectively. 


6^2 + 2j/i + 4^vv^y22/3) 


(A23) 


In Eqs. (A18l-(A23) we have used 


X 3 = 2 — xi — ^2 

1 , 

yi = 2 (1 + y-w - k-x- y-di - X3) 
1, 

y2 = 2 (1 - MlV + Mx “ t^di - X 2 ) 

1, 

ys = ~ ktw ~ y-x + ydi - xi) 


(A24) 
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